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t^^ ■ Abstract: During (re)heating of the Universe after inflation, the relativistic decay 

Q^ ■ products X of the inflaton must lose energy, and additional particles must be pro- 

O . duced, to attain a thermalised state at a temperature Treh- We estimate the rate 

Q ■ of energy loss via elastic and inelastic scattering interactions. Elastic scattering is 

an inefficient energy loss mechanism so inelastic processes, although higher order in 

the coupling a, can be faster because more energy is transfered. The timescale to 

^' produce a particle number density of 0{T^^^^ is the inelastic energy loss timescale 

^ ■ ~ (a^Tj.^jj/i5^)~^. The Universe will be thermalised within a Hubble time at T^ch if 

E^ < a^Mpi. 

X 
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1. Introduction 

The primordial density fluctuations required to account for the observed large-scale 
structure in the universe can plausibly be generated during a quasi De Sitter expan- 
sion phase at early times [1]. The recent detection [2] of the expected character- 
istic signature of such fluctuations in the cosmic microwave background anisotropy 
provides strong evidence in favour of such an inflationary phase during which the 
Universe underwent superluminal expansion while dominated by the vacuum energy 
of a scalar field — the inflaton [3]. This period must last long enough to generate a 
homogeneous universe with small density fluctuations up to the scale of at least the 
present Hubble radius, but must eventually evolve into the hot radiation dominated 
era of the standard Big Bang model. The process by which the vacuum energy is 
converted into relativistic particles can be quite complex. Traditionally only the 
perturbative decay of the inflaton [4] has been considered for studying (re)heating, 
but there may also be non-perturbative transfer of energy to other fields — dubbed 
"preheating" [5]. However, this is not always an efficient energy loss mechanism for 
the inflaton when its decay products have self-interactions [6, 7]. Field theoretical 
studies of the energy transfer from a cosmological scalar field to other particles have 
been performed, but so far only in the context of toy models [8]. Our interest here is 
to determine the time-scale for the thermalisation of the bulk of the vacuum energy, 
when the infiaton can be taken to decay perturbatively into relativistic particles. 

The momentum space distribution of the relativistic decay products (x) of the 
inflaton (0) will not initially be thermal, so these particles must interact with each 
other to redistribute their momenta as well as adjust the particle number density 
to create a thermal distribution. This process is called "thermalisation", and has 
been discussed previously by many authors [15, 16, 17, 18, 19, 20, 21]. There are 
two aspects to thermalisation — reaching kinetic equilibrium, and achieving chemical 
equilibrium. For the first, the momentum must be redistributed among the particles 
present, which can happen via 2—^2 scatterings and annihilations. For the second, 
the comoving particle number density must be modified, e.g. by decays, inverse 
decays, or 2 -^ 3 particle interactions. We assume here that the infiaton decay 
products are initially more energetic and less dense than the distribution to which 
they thermalise. This means that the infiaton decay products must lose energy and 
produce additional particles. This corresponds to generic "new" infiationary models, 
but possibly not to "chaotic" models where infiation occurs at field values beyond 
the Planck scale. 

It has been noted [15, 16] that the infiaton decay products might not be ther- 
malised at Treh, defined as the temperature when they first dominate the energy 
density of the Universe. Thermalisation was estimated in [15] to occur when the XX 
annihilation rate Fann begins to exceed the Hubble expansion rate H, which happens 
well below Trch in many infiation models. A more detailed analysis using the Boltz- 



man equation can be found in [16]. In a numerical study of the thermalisation of a 
gas of semi-classical particles [17] it was found that the particles indeed reach kinetic 
equilibrium after a few hard scattering interactions, i.e. when Fann ~ H . However 
it takes longer to achieve chemical equlibrium since this requires new particle pro- 
duction; it was argued in [17] that the timescale for this is of order a^^ times the 
kinetic equilibration timescale. 

It has been suggested that soft processes can make thermalisation faster [14, 20], 
because the interaction rate for processes with small momentum transfer is larger 
than the hard scattering rate used in ref. [15, 16]. It has also been argued that 
low energy particles can act as a seed for thermalising the energetic inflaton decay 
products, because the cross-section for annihilation with a soft particle is larger than 
with an energetic particle [20]. This scenario of catalysed thermalisation could be 
relevant if particles reach a thermal distribution via 2-^2 interactions and decays. 
Thermalisation via annihilation and decays has been recently discussed, in a Universe 
where the baryon asymmetry is generated via the Affleck-Dine mechanism [21]. This 
is complementary to the present work, where we will concentrate on thermalisation 
via 2 — i> 2 and 2 — > 3 scattering interactions. 

In this paper we discuss which interactions will thermalise a bath of relativistic 
fermions with gauge interactions produced in inflaton decay. The thermalisation 
timescale is usually taken to be an interaction timescale, which begs the question 
"which interaction?". Thermalisation can be very fast, if it happens on the timescale 
of soft scattering processes, because these cross-sections diverge as the momentum 
transfer goes to zero. At the opposite extreme, the timescale can exceed the Hubble 
time, if thermalisation requires hard processes with momentum transfer of order the 
incident particle energy. We argue that the thermalisation timescale is at least as 
long as the timescale for an inflaton decay product, produced with energy of O(m^), 
to lose an energy ~ {m^ — T^eh)- We estimate Feias (= d(lni?)/dt), the rate of energy 
loss via elastic scattering of an energetic particle incident on other energetic particles, 
and find that the timescale for it to lose its incident energy is the hard annihilation 
timescale. We estimate the rate for energy loss due to inelastic 2^3 scattering, and 
find Find ^ Fcias- This suggests that 2^3 scattering interactions can thermalise the 
Universe faster than 2—^2 processes and decays, because although higher order in 
a, 2 — i> 3 interactions are lower order in Treh/""^,^- Such processes are also neccessary 
to bring the particles into chemical equilibrium. The timescale to produce a particle 
number density of 0(Tj.gj^) via 2 — i> 3 interactions is given by T^^cii which is therefore 
the true thermalisation timescale. 

In the next section we introduce our model and review relevant previous work 
[17, 20, 21]. Its purpose is to introduce notation and make this paper self-contained. 
In the following two sections, we discuss what interaction should be used to estimate 
the thermalisation rate. We study the rate of energy loss of a particle via elastic 
scattering in section 3. This process is solvable and only has a logarithmic infrared 



divergence. In section 4, we consider inelastic 2^3 processes. We estimate the 
timescale for an inflaton decay product to lose an energy ~ m<^, and the timescale 
to produce a number density ~ Tj^eh- ^^ section 5, we outline the thermalisation 
discussion of previous papers, which uses the annihilation rate as the thermalisation 
rate. We then repeat the analysis using the inelastic rate, which suggests that the 
Universe will indeed be thermal at Treh- We present our conclusions in section 6. 

2. Model 

We consider a scalar field whose potential energy is the principle component of 
the energy density of the Universe p. (We refer to as the inflaton, but it could 
be any other scalar field, e.g. a modulus, which dominates the universe.) The 
Hubble expansion rate when starts oscillating coherently about the minimum of 

its potential is 

37r 
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The field decays at a rate F^ = a(f,m^ to two light fermions x ^-iid x- ^ The 
energy density of coherent scalar field oscillations redshifts like matter as the Universe 
expands, so the energy density p<^ then decreases with time r as a~^e~^'^. The 
Universe will be dominated by the oscillations until r ~ FT^, when most of the 
energy is transfered to the relativistic decay products. 

The "reheat temperature" T^h is defined when if ~ F^ as 

Prad(arch) = „„ '^^ , (2.2) 

where (?* is the number of relativistic degrees of freedom. The x particles are rela- 
tivistic so we can define a "temperature" T for this radiation, as done above. If the 
X particles have reached kinetic and chemical equilibrium, T will correspond to the 
thermodynamic temperature. 

It is recognised [3] that there is a bath of relativistic particles prior to T^^h since 
decays over time and not instantaneously at Tj-ch- There could be interesting 
implications for baryogenesis and other particle abundances [10, 11, 12, 13] if this 
bath is thermalised. The number density n^ of x and x particles produced as 
decays is 

n^ = 2n^(a,J(l - e-^^^"--)) (^)' , (2.3) 

where n^{ai^) = p^{ai^i)lm^. (This neglects xs pair-produced in x self-interactions.) 
The number density of xs increases rapidly until r ~ 2rin, then decreases as a"^/^ 
until r ~ F7^, when most of the inflaton energy has been transfered to the xs. 



^Our results are largely based on dimensional analysis, so we would not expect tlicni to change 
if the inflaton decay products were different. 



After the xs are produced, their energy redshifts. A x produced at time ti with 
energy m0/2 will have energy E2 = {m^/2){ai/a2) = (m<^/2)(ri/r2)^''^ at time T2, so 
the comoving distribution in energy space is 



AN dA^ dri ^ 2E2TT2 ^, , , , ^1, ,^ ,, 

di^2 dri di^2 \ m^ m^ 

where N{t) = 'n,(r)(a(r)/ain)^. For ti ^ T~^, we can use A'(/,(ti) ~ N^{ai^). The 
energy density in xs at some time r < FT^ will be 

Px(«) = Z'"'^' dE^^ - l^r - nn)pM.) (^)' - ^T^ (2.5) 

The maximum x energy density (which occurs at r ~ 2rin) is defined to be g^n'^T^^-^/ZQ. 
It is easy to see that between Tmax and T^eh Px ~ a"^''^, T ~ a~^/^, and (Tmax/^reh)^ 

We assume that the xs have SU{Nc) gauge interactions among themselves with 
coupling a ~ 1/30. We would like to know how soon the x distribution will have the 
equilibrium form f(k) ~ (e^''^ + 1)~^. We can get a qualitative answer by comparing 
the expansion rate H to interaction rates, in which we make some attempt to include 
factors of tt and N^ in sections 2 and 3. We drop them in sections 4 and 5, where 
the discussion is more approximate. A more accurate result could be obtained by 
solving Boltzman equations, or perhaps other more appropriate equations [9] , for the 
particle phase space distributions. 




Figure 1: x scattering. Time runs from left to right 

There are various 2^2 rates that could be compared to H, such as the anni- 
hilation rate, which is slow, or the scattering rate, which is fast. In the comoving 
rest frame, 2^2 processes redistribute energy because they do not take place in 
the centre-of-mass frame. They can bring a group of particles into kinetic (but not 
chemical) equilibrium. 

The annihilation cross section for xs of energy ~ mfp/2 is 

leA^cft^ . . 

O-ann ^ n ' (2-6) 

ml 



while the cross-section for scattering (see figure 1) is infrared divergent: 
_ ^.iVcVra^ / d(sin^) _ ,_,, o / dt 



2nN,a' / - (2.7) 



ml J {l-cosOy ' J t2 

where 9 is the scattering angle between k and k' and t is the 4-momentuni transfer 
squared. Note that there is no thermal bath whose plasma frequency can provide a 
cutoff for the scattering cross-section (we are trying to compute when the thermal 
bath appears). In the following section, we will show that the thermalisation rate, 
due to scattering, is the rate ~ '^x'^ann associated with annihilations. 

The estimates in this section assume that o;,^ ^ m^/Mpi, and m^ ^ Mgut- The 
assumption that infiaton decay products need to lose energy and produce particles 
to reach a thermal distribution breaks down at a^ ~ m^/Mp\. If we approximate 
the equilibrium number density Ueq ~ T^ ~ p^/^, we find, at T^eh'- 

% ^pT ^ 1 ^4'Mpi ^Z,^^ (2.8) 

so n^ <^ Ueq if a^ <^ m^/Mpi. Perturbative particle interaction rate estimates 
suggest that gauge interactions are not fast enough to be in equilibrium at energy 
densities of C(lO^^GeV)^, so our discussion breaks down at these energies. 



3. energy loss in elastic scattering 

From the perspective of a x particle produced in the decay of a 0, thermalisation 
is a process of losing energy. So the thermalisation timescale for x is the timescale 
over which it can lose energy ~ m^ — Treh ~ ^tT'^p to the surrounding particles. 
We can estimate this timescale by integrating the particle's rate of energy loss via 
scattering. In this section, we consider elastic scattering. We neglect annihilations, 
which we do not expect to qualitatively affect our calculation, because hard scattering 
processes are included. Scattering interactions have been neglected in many previous 
thermalisation papers, who take the 2^2 annihilation rate as the thermalisation 
rate. We show here that the thermalisation rate due to scattering is logarithmically 
enhanced over the annihilation rate. So taking the 2 — i> 2 thermalisation rate Fthcrm ~ 
Tann IS justified, but claims that scattering is irrelevant for thermalisation are not. 
We neglect the cosmological expansion for the following two sections to avoid the 
scale factor cluttering up formulae. 

The elastic scattering cross-section (2.7) diverges as o" oc d6/6^, where 6 is the 
scattering angle. The energy exchanged is AE ^ m^6'^/2, so the energy loss rate has 
a softer infrared divergence. The cross-section for soft glancing scattering may be 
very large, but the energy exchange is small, which reduces the importance of soft 
processes for thermalisation. 



The rate of energy loss of a particle scattering on a thermal bath 

— = (ncTscatAE), (3.1) 

where n ~ T^ is an equilibrium number density and AE is the energy exchanged, 
has been extensively studied in finite temperature field theory [22, 23]. The logarith- 
mically divergent result is 

^r^2a^THog(^^, (3.2) 

dr \ gmin / 

where q^ax and gmin are the maximum and minimum momentum transfer respectively. 
We follow ref. [22] in making a rough estimate of dE/dr for a x scattering on an 
unthermalised bath of inflaton decay products. Consider the kinematics of scattering 
in any frame where p ■ g*^ (or the approximation where p- k' = {k'^/ko) ■ p- k); see 
figure 2 for the momenta of the particles involved. We are interested in the energy 
lost by the relativistic x with incident momentum k. If the 4-momentum exchanged 
is q"^ = (k- k'f = -t, 2 then 

go ^ J ^(1 -cos^) , forp-g = 

ko~\'-^ ,hTp.k' = f^p.k, ^^-^^ 

where (f is the angle between the incident particles, and s = 2kopo{l — cos^p). For 
p' on-shell, go — kot/s. In the relativistic limit where we can neglect particle masses, 
equation (3.1) can be written as [22] 

I dko f d^p f 27rN^a^ (ko - k'o) 

^' ' /(p)(l-cosv5) / dt — — ^, (3.4) 



ko dr J (27r)3''^'' "^V t2 ko 

where we have kept only the most infrared divergent contributions. Here f{p) is the 
momentum space distribution of the xs, and do" = 27r NcO^dt/t^ is the scattering 
cross-section. Using t/s = {ko — kQ)/kQ, one finds [22] 



1 dA;n l67rN,a'^n^ , m^ 27rNca^alM^ 



elas 



"X 



j^m^^ZTTiVeOa^Mpij^m^^ (3.5) 



ko dr ml T^h m^ T^eh 



We took ginin ~ T'reh; we discuss why in the next section. The last equality is 
Feias evaluated at Tj-ch- This is logarithmically enhanced with respect to the hard 
annihilation rate Fann ~ ct^^x/""^^ • 

It is easy to see in equation (3.5) that the infrared divergence of the scattering 
cross-section ~ a^ / dt/t^ is removed from the energy exchange rate because go oc 
t = — g^ (rather than go oc ^/t). Elastic scattering is an inefficient energy exchange 
mechanism, so only speeds up thermalisation by a logarithmic factor with respect to 

'^sA and u arc the kincmatical variables, r is time. 



hard annihilation processes. This agrees with the numerical results of ref. [17], who 
counted the number of scatterings required to bring their bath of particles to thermal 
equilibrium. They found that a few hard collisions per particle were sufficient, but 
that more soft interactions were required. 

We estimate the timescale for a x particle to lose an energy of 0{m^) by elastic 
scattering with other xs to be 

(3.6) 
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4. 2 

There are at least two flaws in the estimate leading to equation (3.6). Firstly, new 
X particles must be created to bring the xs into chemical equilibrium. The average 
comoving energy per x particle will remain m^ until it can be redistributed among 
newly created particles. So the rate of energy exchange, which we calculated, will 
only be the rate of energy loss for an average x if there are newly created x particles 
available to absorb the energy. Secondly, since elastic scattering is an inefficient 
energy transfer process, inelastic processes, although higher order in a, could be 
faster. 

There are various 2^3 interactions by which xs can lose energy and produce 
particles. We do not worry about the spin of the particles created, because this 
section is based on dimensional analysis. We identify particles by their energy — xs 
are inflaton decay products, and gs are the particles of energy ~ Treh being created 
to populate the thermal bath. We refer to gs as gauge bosons, although there will be 
X particles of energy Treh in the thermal bath. We make a separation based on energy 
because cross-sections and number densities depend on it. We neglect gauge bosons 
with energy less than Treh because they would subsequently need to be scattered up 
in energy to attain the equilibrium relation {E) ~ 3T. The x particles can produce 
gauge bosons and lose energy via XX -^ XX9, XX ^ XX 9 and xfi- ^ X9 9- The 
first process is s-channel, with s ~ m^, so we concentrate on the last two t-channel 
processes. 

There are many Feynman diagrams contributing to xX ~^ XX 9-, which can be 
found, with the associated matrix elements, in ref. [24] . We estimate the cross-section 
from the diagram of figure 2 with an outgoing gauge boson of momentum d, {(? = 0) 
emitted from the leg p'. So p^ = {p" + i)"^ = W'^ is off-shell, and p ^ = m^. We 
neglect gamma matrices in the matrix element squared, and estimate 



a,„, ~ 2.iV.„= y -_ ^ log U^ , (4.1) 



dtdW^ a\ (ml 

■ reh \ reh , 

where {k — k'Y = —t, 1/t^ corresponds to the photon propagator, and l/W"^ to the 
fermion propagator. We treat W'^ and t as independent variables, ranging from T^^j^ 



to m'i. We take the lower limit to be T^^^i because we are interested in making 
Ug ~ Tj,^h gciuge bosons of energy ~ Treh- Equation (4.1) might capture the physics 
of the infrared behaviour we are interested in. As t -^ 0, it has the logarithmic 
divergence one expects from bremmstrahlung and a quadratic divergence due to t- 
channel exchange. The remainder of this paper relies on equation (4.1) being a 
reasonable approximation; our arguments will not be true if a^^^^^^^g is suppressed 
by a additional factors of T^c\Jf^<j) due to phase space or cancellations in the matrix 
element. 

One gauge boson for every x can be generated via XX scattering in a time 



Tincl ~ l^yO- 



x^xx^xxgj 



a^T^ 



rch 



m^ 




' elas • 



(4.2) 



To make Ug ~ Tf^^ gauge bosons, we need to make m^/T^eh gauge bosons for each 
X particle, as can be seen from equation (2.8). So the timescale to produce a ther- 
mal distribution of gauge bosons via xX scattering can be estimated as ( jTr^) Tinei- 
However, equation (4.2) is also the timescale for the particles g to thermalise among 
themselves, at a rate ~ a^n^/T"^^^. Once ng > n^, gauge bosons can be produced in 
X9 ~^ X99 and Ug will grow rapidly, as observed in reference [17]. We can write 

drir, 



^ ^ ^«- 



^xxgi^l + ^x^g)- (4-3) 

The first term will dominate until Ug ~ n^ at Tinei, then Ug will grow exponentially 
with a timescale Tinei- So the timescale to produce a thermal number density using 
the cross-section (4.1) is equation (4.2). 

Cooling via the 2^3 scattering cross-section (4.1) will be faster than the 2^2 
cooling rate computed in the previous section. We estimate the energy loss of the x 
which emits a gauge boson to be of order Treh ~ Vt, so 



dE 
d7 



U^X + ^gjO'xX^XXS-'rch/ 



J-' rch/ 



a^{n^ + ng) 



(4.4) 
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The timescale to lose energy ~ m^ through scattering on ^s will be r ~ [a^n^/(m^Trch)]~^- 
This is already a factor oiTj.cii/{am^) shorter than the elastic cooling timescale (3.6). 
However, the xs can lose energy even faster by scattering off the growing bath of 
gauge bosons. If we substitute Ug ~ n^(areh)e'^/'^'"°' (approximately the solution of 
equation (4.3)) into equation (4.4), we find that the timescale for a % to lose energy 
~ m^ is r ~ Tinci- This is the timescale (4.2) to produce a thermal bath, and equals 
^reh/('^"^0) ^ ''"elas- We expect aiJi^/T^^Yi > 1, so 2 ^ 3 interactions cool the xs faster 
than elastic scattering, despite being higher order in a. 

There are two disturbing features to our estimate of the thermalisation rate 
Tinei = Tjnei- It is larger than the elastic rate, which is lower order in a, and it is 
infra-red divergent- We can estimate the cooling timescale due to 2 ^ 3 processes 



of an energetic particle incident on a thermal bath. In this case, dE/dr due to 
2^3 scattering is ~ a^T"^, an 0{a) correction to equation (3.2). This is reassuring, 
because this thermal energy loss rate has been carefully studied [22, 23]. We find a 
faster thermalisation rate for our inflaton decay products at 0{a^) than at 0{a'^) 
because the energy exchanged and target number densities are larger in the inelastic 
case. Our inelastic rate is one factor of m^/Ty-ch larger than the elastic rate because 
Hg > n^; the x scatters more frequently off the bath of particles it created in earlier 
interactions. It is a second factor of m^/T^ch larger because more energy is transfered 
in an inelastic collision than in an elastic one. We have two small parameters in our 
reheating problem: a and Ty-eh/tri^. The inelastic thermalisation rate is higher order 
in a than the elastic rate, but lower order in T,.f.h/^n<f>- 

Equation (4.4) has an infra-red divergence: the rate for cooling by emission 
of gauge bosons of energy ^ diverges as l//i^. Physical observables should not be 
infrared divergent, so a cutoff is required for our estimate. There is initially no 
thermal bath present to justify using the "thermal mass ~ gT" — but there will be 
at the end of (re)heating, so we imagine that the final state gauge boson must have 
energy of order Treh to be on-shell. To see this, suppose that in a time Ar a fraction 
/x^/n^ of the xs scatter inelastically, emitting a gauge boson of energy /i. This creates 
a bath of gauge bosons with number density ng{fi) ~ /i^, through which the next 
generation gauge bosons must propagate. So these next generation gauge bosons 
must have energy ^ fi. As the time interval Ar lengthens, the cutoff /i grows to 
Treh- We therefore require the final state gauge bosons to have energy ^ Treh, and 
assume that the momentum transfer in the scattering is of order the energy of the 
emitted gauge boson. The exchanged gauge boson may not have a "thermal mass" 
due to interactions with other gauge bosons, because it may not live long enough to 
interact with them. 



5. Thermalisation 

We would like to know the thermalisation timescale Ttherm after reheating. We esti- 
mate Ttherm = rtherm) ^.ud identify the thermalisation rate as the rate of energy loss 
for a Xj ot the rate of particle production. (The two turn out to be comparable). In 
this section, we estimate when the Universe will be thermal, using the rates from the 
previous two sections. In the first part, we discuss two-to-to rates [15, 16, 20, 21]. We 
review the thermalisation mechanism suggested in reference [20], and explain where 
we disagree with those thermalisation estimates. We make some remarks on using 
decay processes as a particle production mechanism. In the second part, we compute 
the upper bound on the inflaton mass, below which the Universe will be thermalised 
at Treh due to 2 — > 3 scattering processes. 



5.1 XX -^ XX 

The elastic rate (3.5) (without the log term , which we also drop) has been taken as 
the thermalisation rate in some previous work [15, 21]. This assumes that a thermal 
number density is rapidly produced by the decays of particles involved in the 2 — i> 2 
processes. The rate scales after Treh as Oreh/o- We start in the instantaneous decay 
approximation to review thermalisation estimates based on equation (3.5) [15, 20, 21]. 
Thermalisation will happen immediately at Tj-eh if 

2 

For m^ of (9(10^) GeV, the bound (5.1) translates into a^ > 10"'^*' which is certainly 
satisfied. However, for m^ ~ Mgut ~ lO'^^GeV, a^ > 10~^. For an inflaton mass of 
order the hidden sector scale lO^^GeV, equation (5.1) implies that the inflaton decay 
products can annihilate with each other within a Hubble time if a^ > 10~^^, which 
may not necessarily be the case [15, 21]. For instance, if decays gravitationally, 
a^f, ~ m^/Mpi ~ 10^"*^^. If the (pxx interaction is of electron yukawa strength, then 
a^, ~ 10-^1 

Equation (5.1) appears peculiar, because it gives a lower bound on a^. At 
smaller a^, the inflaton decays later, so the number density of inflatons is redshifted 
and therefore the number density of xs is similarly smaller {n^ ~ n^). The energy of 
the xs remains m^/2, so the x interaction rate is smaller (scales as a^, see equation 
(3.5)). The expansion rate H is also smaller when a^ is smaller, but this is a less 
important effect because H ~ a^m^ at Treh- 

We now relax the instantaneous decay approximation and consider whether the 
Xs produced between T^ax and Treh have time to annihilate. Repeating the estimate 
that lead to equation (3.5), using n^ from equation (2.3) with (1 — e^^"^) ~ Fr, and 
taking E^ = m^/2 gives 

F ,,^^, 32iVea^3iJ^(a)Mgl 2F^ 
'^^^^^^- ml 8nm^ 3H{ay ^^'^^ 

If equation (5.1) is satisfied, Feias > H from Tmax onwards; that is, if Feias > H is 
true at Tmax if it is true at Treh, and vice- versa. Fcias/-f^ is a constant between Tmax 



rsj 



a 



-3 



and Treh because Feias scales as n^ ~ a~^/^, and H^ scales as p^ 

Now let us suppose that equation (5.1) is not satisfied. As discussed in refernce 
[20] , the annihilation rate among the first xs produced at Tmax grows relative to the 
expansion rate H. So it is claimed [20] that these particles can interact and produce 
particles at some intermediate scale factor Cmax < a < a^^h- More energetic particles 
can then thermalise rapidly, when they are produced, by interacting with this soft 
tail after it has reached kinetic and chemical equilibrium. This scenario is dubbed 
catalysed thermalisation [20]. However we disagree with some of the estimates in 
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reference [20] , because the annihilation rate of a x of energy E2 with less energetic 
Xs is always smaller than the annihilation rate with more energetic xs- We can see 
this by evaluating those two rates at T^ch'- 

dE—-——< dE— -— — ~48A^c TrnJarch)- (5.3) 

dEEE2 Je2 dEEE2 m^E2 ^^ ' ^ ' 

We agree that the less energetic xs are more likely to interact (the right hand side 
of equation (5.3) is larger than that of equation (3.5)), but they will annihilate with 
one of the energetic xs whose density is higher, hence their energy will increase. It 
will not substantially cool the energetic xSj because few of them interact with a less 
energetic x- Nonetheless it is possible that a thermalised seed is formed out of the 
less energetic xs. Thermalisation is a process of energy loss and particle production, 
which are both assumed to occur rapidly due to decays in the annihilation and decay 
scenario of [20, 21]. Suppose that x cannot decay, but can annihilate into some 
rapidly decaying particle p. An efficient cascade decay of p could then produce 
many lower energy particles that rapidly thermalise. The mass of at least one of 
the particles at every step of the cascade must be heavy enough that their decay 
rate F^ = X^m'^/E^ is larger than H. This requires some tuning of the mass, and a 
sufficiently large A, but is certainly possible, particularity in supersymmetric theories 
where ffat direction vevs could provide such an intermediate scale mass. We do not 
discuss further the scenario of thermalisation via annihilations and decays [20, 21]. 
In the next subsection, we focus on the generation of particles by 2 ^ 3 gauge 
scattering interactions, which should occur generically in all models. 

5.2 2—^3 interactions 

In section 4, we estimated the energy lost by a x scattering and emitting a gauge 
boson to be Finei ~ Q;(m0/Ti.ch)^rcias. One factor of (m^/T^ch) arises because we 
took the inelastic energy loss in a collision to be Treh, rather than T^^yi/''^<P ^^ ^^ the 
elastic case. The second factor of m^/Tj-ch comes from scattering the xs off the denser 
thermal bath being produced in the 2^3 interactions. We expect that a ~ 1/30, 
and that thermalisation would take place within a Hubble time if T^^h ~ ^^0/6. So 
Q;m^/Tj.gj^ > 1, energy loss via 2 — > 3 processes is faster than via 2^2 elastic 
scattering, and we estimate the thermalisation rate at Treh to be Ftherm = Finci = t{^Ii 
(see equation 4.2). This coincides with the rate to create a number density Ug ~ T^^^ 
of gauge bosons via XX ^' XX 9 and XQ -^ X9 9- Thus Finei will exceed the expansion 
rate at T^eh if 

m^ < a^Mpi. (5.4) 

The COBE results require the inffationary scale (and hence m^) to be much less 
than 10^^ GeV [1], so this condition ought to be satisffed. As expected, this is a 
weaker bound than equation (5.1). It is independent of a^ because both FincKorch) 
and H{aT-ch) are proportional to o;,^. 
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Note that 2-^3 interactions can thermalise the relativistic particles at T^ax 
within a Hubble time if a^ > 'm|/(Q;^Mpi), which is similar to the condition (5.1). 

The rate of energy loss of an energetic particle incident on a thermal bath is 
dE/dt ~ a^T^ (see equation 3.2) [22, 23]. The thermalisation timescale for a particle 
of energy m^ is therefore ~ (a^T^/m^)"^. Our estimate of 2 — > 3 interaction rates 
suggests that reheating a cold Universe is a factor of a slower; the inflaton decay 
products rapidly produce a bath of soft particles, and cool by interacting with them. 

6. conclusion 

If the inflaton decays perturbatively in a cold Universe, its decay products must 
interact to thermalise. There are two aspects to thermalisation: producing additional 
particles and distributing energy among them, so as to obtain a kinetic and chemical 
equilibrium distribution in momentum space. It is often assumed that thermalisation 
proceeds by annihilations and decays, so the thermalisation timescale is taken to be 
the timescale of hard annihilations among inflaton decay products. In this paper 
we discussed thermalisation via scattering interactions. Our estimates suggest that 
these soft processes lead to faster thermalisation. 

We first considered energy loss via scattering between inflaton decay products 
X, and found that the timescale for a % to cool down to Tj-eh is the timescale of hard 
annihilations. The scattering cross-section is infrared divergent so the interaction 
rate is large, but little energy is exchanged in soft elastic scattering so it does not 
lead to rapid thermalisation. 

We then estimated the cooling rate via inelastic 2^3 processes Finei, and found 
it to be much larger: Finei ~ am^/T^gj^x the elastic rate. We also estimated the 
timescale to create a particle number density ~ T^^^ via 2-^3 interactions and found 
that it was of order T^^^. Thus 2^3 scattering is a more efficient thermalisation 
process than 2^2 elastic scattering among inflaton decay products, because the 
energy transfer in a collision is a factor of m^/Treh larger, and because inflaton decay 
products can collide with the particles they produced in earlier 2 — ;> 3 interactions. 

Our estimated thermalisation timescale at T^eh is 






Tthcn-r^ - P£^ (6.1) 



where n^ is the inflaton number density just before Tj-ch- This result is independent 
of the details of the reheating model. If the Universe is dominated by relativistic 
particles of energy E, then they will thermalise within a Hubble time to a temperature 

Treh ~ P'/' if 

3 

Ttherm -^ ^^ > H ^ E < a^ Mp^. (6.2) 
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